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(«~^ , Abstract. We present preliminary results of study of the r|7t ,r\ n and ^7i7t-systems produced in 

■ the Tt^Be-interaction at 28 GeV /c. j''^ = 2++ and l^+-waves resulted from the PWA have been 

psj ' fitted in each system separately to establish the nature of the 1 ^-wave. A hypothesis of the 1 ^- 

^ . wave resonant nature in the r|7t^ and r|'7i^ has no statistically significant preference over the non- 

O ■ resonant one. The ^JiTt-system analysis confirms the results of the 37 GeV / c-heam data analysis in 

I favor of the resonant treatment of the bump at 1 .6GeV. 

o : 

cn . 

INTRODUCTION 

>: 

lO ■ At present several groups have evidence for 1 ^ meson production, which is forbidden 

^ ! for ordinary quarkonia and hence might be a good candidate for a hybrid [@]. The 

I 7ri(1400)-meson shows up in the r\n final state in two experiments VES 

^ ■ group found that the results of the r\'n^ and Z^iTT-system Partial-Wave Analysis(PWA) 

^ ! at 31 GeV /c agree with the production of the higher mass 7ri(1600) evidence for 

which in the r\'n^ has been also confirmed by E852 at BNL [0]. A broad bump observed 
by VES in the l^+-wave intensities of the r['n, biK, pjl at l.6GeV was assumed to be 
a single object. The resonant nature results from the analysis of the J^^ = 2++ and 
l"+-waves in the qoti^'kP [§]. 
^ ' A data sample was collected in the 1996 data run at VES spectrometer exposed 

> . by the 2SGeV/c momentum TT^-meson beam. A detailed description of setup can be 

^ i found elsewhere [|^]. Over 5x10^ triggers were recorded during the data taking period, 

' which is ~ 2.5 more than at 31GeV /c. A trigger is an event with at least two charged 

tracks from the beam interaction moving in the beam direction and falling into the 
spectrometer aperture. Events of the studied final states production form subsets of 
reaction n Be — > 7l+27l^ +^7r^(ri) +Be, where k= 1,2. To select events of this reaction 
the following criteria were applied to the reconstructed events: 

• 3 reconstructed tracks should form the TT+TT^TT^-system; 

• an interaction vertex must be inside the target; 

• we require for 2-4 y-clusters in the calorimeter; 
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• a 2Y-system with mass closest to the mean 7r (ri) -meson mass within 30(78)Mey is 
identified as the 7r^(ri) -meson. A IC-fit to the Tl'^{r\) mass is applied to the identified 
YY pairs; 

• total visible energy Efot must satisfy the following in-equation: 25GeV < Efot < 
30GeV. 



STUDY OF THE riTi" -SYSTEM 



A decay chain r\ n'^n~Tr, tt — > 2y was chosen for the titt" production study. A 
clear peak of the Tj-meson corresponding to the riTC production is observed in an 
invariant mass spectrum of the TT+TT^TT^-system in fig. [l]a. The shape of the peak is 
parametrized with a sum of two Gaussian functions and resolution of the narrow one 
is equal to 6MeV. Shown in fig. [l]b is the invariant mass distribution of the n^ln^'KP 
for events from the r\ region (530 < M^+^-^q < 566MeV) with superimposed plot 
for sidebands (502 < M^^+^^-^o < 520Mey and 576 < M, 



< 59AMeV) which are 
used for the background estimation. The background subtracted riTT invariant mass 
spectrum (fig. [I]c) shows a dominant (32 (1320) -meson production with a small bump 
near M^^^j- 



IGeV of a possible (980) -meson production. 
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FIGURE 1. Effective mass of the ti+ti ti" a), 7I+27I %^ b) and riTi c). — f ' distribution d). 

A negative value of the invariant four-momentum-transfer between the beam and the 
riTT^ squared — ^'-distribution in fig. |l]d shows the characteristic features of the projective 



wave production. Here t' = t — tmin, where is a minimum transfer squared for a given 
TjTi;^ mass. 



The T] 71 -system PWA 

We applied the PWA procedure to expand the data in terms of the partial waves. An 
event is considered in the spirit of the isobar model [ 1T0| ] as a production process followed 
by a chain of the subsequent decays into the ri7i~ and r\ tt+tt^tt^. The production pro- 
cess of the riTT^ is assumed to be described by a rank one density matrix for each natural- 
ity [ pT] ] . An event probability for each naturality is equal to an average over possible non- 
interfering Tt^nYvP combinations Ncomb of a product: Pg^, = jj^Y^iPiDiim-iTAGiimj,-^) 

^ ^^comb 

due to presence of two tt^'s, where: - is a probability of the T] — > Tt^njn^ decay 
which corresponds to the decay Dalitz-plot distribution, - is the Tj-meson shape pa- 
rameterization being used to take into account the Tt^Tt^n^ mass resolution. Here P, is an 
expansion squared of the product of the two pseudo-scalar production and decay ampli- 
tudes in terms of the partial wave decay amplitudes. The decay amplitudes are defined 
by the following set of quantum numbers J^M^, where J stands for a total angular mo- 
mentum and is equal to an orbital momentum L in the system of two pseudo-scalars. 
With P we denote parity, M - an absolute value of the projection. Here r\ is exchange 



naturality. The decay amplitude is parametrized with spherical harmonics [ |12| ] multi- 
plied by a breakup momentum involved in the L-th degree. Under the assumption about 
the rank of the density matrix the production amplitudes are denoted as: Lq for the M = 
and for the M = 1 waves. 

The complex production amplitudes are determined from the extended maximum 
likelihood fit [|13p. 

A wave set includes 3 subsets of partial waves non-interfering with each other. The 
first is composed of waves with negative naturality, the second - with positive one and 
a subset made of a single wave FLAT which describes the background in the 7r+27l~7l*^- 
system. 

The ambiguous solutions [ [T^ ] were found by repeating 100 times the fit from the 
random starting values of the fit parameters. 

50300 events in the range 0.8 < M{n+2n^n^) < lAGeV, \t'\ < IGeV^, with the mass 
of at least a single n^n^n^ combination in the T] region or the sidebands were subjected 
to the PWA. The PWA has been carried out independently in each of the AOMeV bins. 
To make the data set cleaner the following additional cuts were applied: cos0^„' < 0.8 

and the track angular separation in projections AAX{AY) > 10{6)mrad, where &'^o is 

the "kP polar angle in the overall CM helicity reference frame. 

In fig. are shown predicted by the PWA fit acceptance-corrected wave intensities 
and a relative phase between the P+ and waves as a function of the riTT^ mass. With 
a thick line is shown a range of possible ambiguous solutions and with a thin one the 
maximum extent of errors. 

The intensities of the S and P waves with the negative naturality are small and 
consistent with zero in all the M^^^- range. The Dq wave has a non-zero contribution 
in the M^^^- > 2GeV. We found a significant D intensity in the whole mass range of 




unclear nature. 

The D_|_-wave is the dominantly produced wave in the V['K~, where a huge peak 
corresponding to the 02(1 320) production is observed. A broad bump with a spike at 
M^j^- = lAGeV is observed in the P+. The spike might results from the fit quality 
worsening in the a2(1320) region and where one should notice the increase of errors. 
A sign ambiguous relative phase shows a rapid motion in the 02(1320) region. 

We have tried also to include the waves with the orbital momentum higher than 2 and 
the waves with M > I, but they were found to be not significant. The quality of the fit 
is controlled by comparison of distributions of variables which describe the T|Ji~ -system 
for the experimental and Monte-Carlo data. The Monte-Carlo events were generated 
according to the described above model with the values of the production amplitudes 
resulted from the best fit solution in a mass bin. 



The TjTi -system mass-dependent fit 



In order to establish the nature of the P^-wave we have carried out a combined mass- 
dependent(MD) fit of the expected production amplitudes to the P+ and Z)+ resulted 
from the PWA in the whole mass range. Fitted with the method are the and 
intensities and the real and imaginary parts of the production of both amplitudes with 
the corresponding errors in each mass bin. 

We assume two general production manners of a resonant and a non-resonant type. 
The resonant amplitude is parametrized with a Breit-Wigner form: 



where: m - the riTT mass, mg - a resonance mass, Fq - a nominal width, F^ - a partial 
width in the final channel, Ttot - a total width, q{qo) - a breakup momentum of the r\'K 
with the m(mo) mass, - a barrier factor [|T3|], BRf - a branching ratio into the f-th 

final state. 

By the notion background we just mean the non-resonant production manner with- 
out going into the nature of the process. The main difference of the background 
type from the resonant one is absence of phase motion. The background amplitude is 
parametrized with the following free form Ag — LIPSi x (m — m;)"exp (— P(m — nit)), 
where: LIPSi{m) - a phase space, - a threshold mass, a and |3 are shape parameters. 

A procedure of the wave construction from the objects is in the following. At first 
a wave is formed of evident structures such as the a2(1320) in the D+-wave. If there 
are no such objects we take an object of any production type arbitrarily. After doing the 
fit to thus constructed wave we then necessarily try to describe the data with the object 
of the alternative type. At the next steps we add consistently another objects of both 
types to both waves until we will succeed in the description of the wave intensities and 
their interference in the whole mass region. All objects are added coherently with their 
own mass -independent phase. The parameters of the MD fit are: the mass, width and the 
number of the resonance events, the shape parameters and the number of the background 
events, the phase relative to a reference object. 



TABLE 1. Compai-ison of the MD fit vari- 
ants 
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To begin with we have fitted the a2(1320) Breit-Wigner to the intensity and failed to 
describe the wave at masses > l.SGeV. Then we have applied the described above 
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procedure. At first the was fitted to a sum of the a2(1320) and the background and 
the to the background. The second fit was carried out with the wave background 
replaced by the resonance, which will be called further as an aj-meson for convenience. 
There have just been reported the evidence for the (32(1750) [|T^. The next steps were 
to fit the wave to the resonance, called a 7ri(1400). The evidence for it have been 
reported by E852 and CBAR groups [|]], 0], [|]] and which production did not contradict 
to VES results [[T7|]. The results of the MD fits in the case, when the has been fitted 
to a sum of the ^2(1320) and the background and the P^ to the background or to the 7li 
are shown in fig. ^ with thick solid and dashed curves respectively. 
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FIGURE 3. The r|7t D+ a) and P+ b) intensities. The real c) and imaginary d) parts of D+P+ product. 
The D+/P+ relative phase e). The smooth curves are the MD fit results. 

The summary of the MD fit results are presented in tab. [I]. The objects which make 
up the fitted partial waves are marked by the plus sign in a row. A %^ over the number of 
degrees of freedom of the corresponding MD fit is shown in the fifth column. 

From the comparison of the first four rows one can conclude that there are no statis- 
tically significant difference between a single resonant and non-resonant object descrip- 
tion of the P+. 

The fits of the P+ to a sum of the TTi and the background are shown in the last two rows. 



A fit with free TTi parameters resulted in the following values: M, 



a2(1320) 



m6±\MeV, 



ra2(i320) = 1 13 ± 3MeV, and M^.^ = 1316 ± 12Mey, T^t, = 287 ± 25MeV. The quoted 
errors are statistical only. Statistical significance is nearly the same as for the fit with the 
TTi parameters fixed to the values reported by E852 A closeness of the TTi mass to the 
02(1320) mass may indicate that there is an underestimated influence of the 6(2 (1320) 
on the 



STUDY OF THE ti'ti" -SYSTEM 



A decay chain r\' r\K^n^, r\ 2y was chosen for the ti'ti:^ production study. A clear 
peak of the r\' corresponding to the ti'tt production is observed in the invariant mass 
spectrum of the r\n^n^ in fig. Shown in fig. ^ is the invariant mass distribution of 
the r\n^2n^ for events from r\' region (940 < M^^^+t^- < 916MeV) with superimposed 
plot for the ri' sidebands (912 < M^^^+t,- < 930MeV and 986 < Mi^jt+ir- < lOOAMeV) 
which are used for the background estimation. The background subtracted Tj'TT invari- 
ant mass spectrum (fig. |^d) shows a clear a2(1320) peak and a broad large bump at 
Mf^j^- = \ .6GeV, which strongly differs from the corresponding distribution of the riTT". 

X 10 




FIGURE 4. Effective mass of the yy a), ri7t+7t b), ri7t+27t c) and ti'ti d). — f ' distribution e). 
A —t'^i^- distribution in fig. ^ shows characteristic features of the projective wave 
production like for the riTT . 

10700 events in the range 1.15 < M(ri7r+27r") < l.\5GeV, \t'\ < IGeV^, with the 
mass of at least a single rjTr+Tr^ combination in the Tj'-meson region or the sidebands 
were subjected to the PWA. The PWA has been carried out independently in each of the 
50MeV wide bins. 

We have applied the same procedure as for the riTT^ -system PWA, described in the 
appropriate section. In fig. ^ are shown predicted by the PWA fit acceptance-corrected 
wave intensities and a phase difference of the P+ and D+ as a function of the r\'Tl^ mass. 
With a thick line is shown a range of possible ambiguous solutions and with a thin one 
the maximum extent of errors. 

The intensities of the negative naturality waves are small and nearly consistent with 
zero in all the M^^/jj- range. The Dq wave only shows some non-negligible contribution 
in the regions lA^1.6GeV and > 1.9Gey. 

The 02(1320) peak is observed in the intensity. But contrary to the riTl^-system 
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FIGURE 5. The r|'7i wave intensities and the D+/P+ relative phase. 



behavior this wave shows a broad structure in the region M^i^^- > l.SGeV. The rate of 
this high mass production relative to the a2(1320) is larger for the r\'Ti^ than for riTT^. 

A broad bump of a ~ 350MeV width is observed in the at M^f.^- = l.6GeV. A 
sign ambiguous relative phase shows a rapid motion in the 6(2(1320) mass region. 

We have tried also to include the waves with the orbital momentum higher than 2 but 
they were found to be not significant. The quality of the fit is controlled by comparison 
of distributions of variables which describe the ri'TT^-system for the experimental and 
Monte-Carlo data. 

To establish the nature of the P+-wave we have carried out the MD fit following the 
procedure which had been applied for the riTT^ study. The partial wave amplitude 
was parametrized with a sum of the a2(1320) Breit-Wigner function and the background 
or an This additional object was included to describe the broad bump in the region 
^r\'Tz- > 1 -SGeV . The partial wave was parametrized with the background or with the 
711(1600). The indication to the possible resonance nature of the bump at \.6GeV was 
reported by our group as a result of the r\'n^ PWA at 31 GeV /c . Later the observation 
of the resonance have been reported by E852 [0]. The result of the MD fit in the case, 
when the is fitted to a sum of the 02(1320) and the a2 and the P^ to the resonance 




only is shown in fig. ^ 

We have found no statistically significant preference of the fit with the P+ saturated 
with the resonant over the non-resonant description. In the case of the partial wave 
amplitude parameterization with a sum of the resonance and the background the statis- 
tical significance is larger than in the previous cases. However the a2(1320) width was 
found to be somewhat larger than it could be expected from the resolution study. In the 
cases, when the is composed of the a2(1320) and a a'2, while the includes the n\, 
the a'2 mass was found to be somewhat 1.8 ^ \ .9GeV and the width OA^O.SGeV, while 
one might expect to find the same parameters as for the a2(1750). 

The above consideration of results prevents us from drawing the unambiguous con- 
clusion about the nature of the P^-wave in the ti'tt^. 



STUDY OF THE COTI TI^-SYSTEM 

A decay mode co n^n^n^ was chosen for the oon^'nP production study. A clear peak 
of the CO meson corresponding to the Qon^n'^ production is observed in the invariant mass 
spectrum of the 'K'^'K~'nP in fig. Shown in fig. ^ is the invariant mass distribution of 
the n~^2-K~2-nP for events from the (O region (758 < M^+^-^o < SOSMeV) with super- 
imposed plot for the (O sidebands (684 < M^+^-^^u < l\9MeV and 846 < M^+^-^a < 
SllMeV) which are used for the background estimation. The background subtracted 
COTT^TT^ invariant mass spectrum (fig. ^) shows an a2(1320) shoulder and a broad large 
bump atM^j^-jjO = l.SGeV. 




FIGURE 7. Effective mass of the 7t+7t Tt" a), 7t+27i 271" b) and (an Ti" c). -t' distribution d). Effective 
mass of the cop" e) and Tt^Ti" f) from the coTi^Tr" for 0.08 < |f'| < Q.lGeV^. 



A —?^j^_j^o -distribution in fig. ^ shows a sharp peak in the low \t'\ < O.OSGeV^ (LT) 
region, corresponding to the coherent production on a nucleus, which flattens in the 
high \t'\ region (0.08 < \t'\ < OJGeV^) and which wiU be called further the HT. This 
behavior in the HT-region corresponds to the incoherent production on a nucleon. The 
projective wave production dominates in the HT-region. Therefore we have divided a 
data sample into the LT- and HT-sets, which were analyzed separately. A large p^- 
meson peak is observed in the invariant mass spectrum of the ti^'kP, produced with the 
CO(see fig. ^jf). This means the dominant production of the (Op . There is a peak of the 
Z7i(1235)^-meson in the oon invariant mass spectrum shown in fig. corresponding 
to the bi{l235)n production. The same is for the l?i(1235)°. 

A PWA of 284000 events in the range 1.2 < M(7r+27r"27r°) < 3GeV, 0.08 < \t'\ < 
OJGeV^, maxjj±(cos0'^'^') < 0.92 with the mass of at least one 7r+7r^7r*^-combination in 
the (0 region or in the sidebands has been carried out. Here 0''^' is a track polar angle in 
the overall helicity CM reference frame. Events in each of the 50MeV or lOOMeV wide 
bins were fitted separately. The Illinois method [ [Tl| ] was used in the PWA, however with 
an important improvement. We added the possibility to restrict a rank of the density 
matrix [|18]]. The partial waves are denoted by the J^M^LS {isobar — bachelor), where 
S stands for a total spin in the isobar-bachelor system. The partial waves of the (Dp , 
Z7i(1235)7r, p3(1690)7r, pi(1450)7r are included in the wave set. The results of the PWA 
are consistent with our previous analysis of the cok^'kP at 31GeV /c [|T9p. A 2+l+52(a)p) 
shown in fig. ^ was found to be a dominant wave with a clear a2(1320) peak and a 
broad bump at M^^-^o = l.lGeV . A significant l^l^ S\{b\n) wave shown in fig. |B|b is 
observed with a broad bump at M^^^-^ = l.6GeV, which is in the maximum ~ 15% 
of the ^2(1320) peak height. The important feature of the combined 2+ and 1^ wave 



behavior demonstrated in fig. |8|c is a region bounded to 1.5 < M^^-^o < 2GeV, where 
the coherence is significantly non-zero. It should be noticed in fig. ^ an 80° rise of 
the l~-wave phase relative to the 2+-wave phase right in this region and which may be 
attributed to a 1 resonance. 

The MD fit has been carried out to establish the origin of the l~l"'"-wave in the bin. 
The high intensity 2+l+52((Op) was selected to be a reference wave. The fit has been 
done in the same way like the r\n^ fit. However absence of unambiguous knowledge 
about the nature of the bump at M^^-^o = 1 JGeV and an extension of the density matrix 
rank to three in the PWA requires inclusion of additional fitting parameters in the MD fit. 
The naive model is to describe the 2+ wave by a sum of the a2(1320) Breit-Wigner and 
something, which has a form of the bump at IJGeV, while the 1^ may be constructed 
from a broad incoherent background and from a resonance, which is coherent to the 2+ 
bump. 




3000 - 




2000 



1000 



2 : 




e) 








-100 - 



FIGURE 8. The 2+l+52(cop) a) and 1^1+51(^171) b) intensities. A coherence parameter c) [HI]. The 
real d) and imaginary e) parts of their non-diagonal p-matrix element. The 1^ phase relative to 2+f). The 
smooth curves are the MD fit results. 

The fit describes the interference pattern satisfactorily as seen in fig. |[ We have 
tried another hypotheses: to fit the 1^ wave to a partially coherent background or to 
make up the 2+ wave from a a'j in addition to the a2(1320). However we failed to find 
qualitatively better solution. The data are consistent with the resonant description of the 
l"l+51(Z?i7r) with the mass 1.6GeV and the width 0.33GeV. 



CONCLUSIONS 



We have carried out the PWA of the riTT ,r\'n , Z?i7r systems produced in tt 5e interac- 
tion at 28 GeV/c. 



We have tried to understand the nature of the observed partial waves in each system 
by the fit to the resonant and non-resonant description. 



We observed the P+ in the r\n^ which exhibits the phase motion with respect to the 
in the coherent model. We could not unambiguously establish the nature of the 
P+ in the riTT^. 

The production of the ti'tt is dominated by the P+ and D+. The common analysis 
of the observed Tj'TT^ partial waves gives no preference to the hypothesis of the 
7li(1600) production. 

We confirmed the results of the bin analysis of the SYGeV/c-beam data. The broad 
bump has been observed in the l^l^Sl{bin) of the a)7i~7i'^. The simultaneous fit 
with the 2+l+52((jop) favors the existence of a resonance in the bin with the mass 
~ 1.6GeV and the width ~ 330MeV. 

The PWA of the 31GeV / c-beam data sample in the several t' intervals shows a wide 
bump in the J^M^ {isobar - bachelor) = l"l + (p7T:)-wave atMp^ < IGeV [Q. 
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